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Impurity effect is systematically studied in doped Haldane material Pb(Nii_3;M2:)2V208 (M = 
Mn, Co, Cu, and Mg) by use of DC and AC susceptibility, and heat capacity measurements. The 
occurrence of three-dimensional ordered state is universally observed for all the impurities and the 
complete temperature - concentration phase diagrams are obtained, which are qualitatively similar 
to that in other spin-gap materials. The unique feature is found in the drastic dependence of the 
transition temperatures on the species of the impurities. The consideration of effective Hamiltonian 
based on VBS model makes it clear that the ferromagnetic next-nearest-neighbor interaction and 
the antiferromagnetic nearest-neighbor interaction between impurity and edge spins play a key role 
in the unique feature. 



I. INTRODUCTION 

Low-dimensional quantum magnetism is one of the 
most exciting topics in condensed matter physics. The 
simplest model is realized in one-dimensional Heisenberg 
antiferromagnet, where the Hamiltonian is expressed as 



J > 0. 



(1) 



Exact ground state in S' = 1/2 is non- magnetic singleti 
and its excitation is gapless known as des Cloizeaux- 
Pearson mode.^ In 1983 Haldane pointed out a quali- 
tative difference between spin integer and half integer 
in Eq. JQ); in the former spin correlation decays expo- 
nentially and spin excitation has a spin gap, while in the 
latter the correlation decays by power law and the excita- 
tion is gapless^i^ A mimic model of spin integer Heisen- 
berg antiferromagnetic chain known as AKLT (Affleck, 
Kennedy, Lieb, and Tasaki) Hamiltonian^'^ was proposed 
by adding a quadratic term in Eq. AKLT Hamil- 
tonian can be solved analytically and the ground state, 
which is known as valence-bond-solid (VBS) state, is sim- 
ilar to that of simple Heisenberg Hamiltonian. In a finite 
chain staggered spins are induced in the vicinity of its 
edges^ and they can be approximated as a couple of 
effective S — 1/2 spins on edges. Edge spins predicted 
by AKLT Hamiltonian is realized in the impurity-doped 
Haldane chain, as is schematized in Fig. ^(a), and the 
low energy excitation can be expressed by the effective 
Hamiltonian consisting of three spins with open bound- 
ary condition; 



H = Jm(si-S_ 



M 



S2). 



(2) 



The predicted energy excitations were observed in the 
ESR spectrum in impurity-doped Haldane material 
NENP (Ref. 0) and numerical calculation also supports 
application of AKLT Hamiltonian to Haldane chains 

PbNi2V208 is a new member of inorganic Hal- 
dane materials discovered by Uchiyama et al. in 
1999»i2i The Haldane gap was confirmed by the 



(a) Haldane chain + impurity 



^1. 




Ni''(S=1) 



Imi 




purity M 



(b) Haldane chain + NNN interaction + 
non-magnetic impurity 




(c) Haldane chain + NNN interaction + 
magnetic impurity 




FIG. 1: (a) Effective spin configuration in doped-Haldane 
chain. 5 = 1 spins are described by a couple of 5* = 1/2 
spins (filled circles) , and symmetrization operator (ellipsoids) . 
Adjacent 5 = 1/2 spins on different Ni^'^ sites are coupled 
by singlet bond (solid line), (b) Effective spin configuration 
in non-magnetic impurity doped Haldane chain with next- 
nearest-neighbor interaction. The next-nearest-neighbor in- 
teraction Jnnn and nearest-neighbor interaction Jnn is indi- 
cated by thin and bold solid lines, respectively, (c) Effective 
spin configuration in magnetic impurity doped Haldane chain 
with next-nearest-neighbor interaction. 



bulk susceptibility, high-field magnetization, and in- 
elastic neutron scattering. ^"■^"'^ The unique feature in 
PbNi2V20g is spin- vacancy-induced antiferromagnetic 
phase in Pb(Nii_j;Mgj;)2V208r2*i^ where coexistence of 
gapless and gap excitations is suggested. The occurrence 
of the antiferromagnetic phase is attributed to the fer- 
romagnetic next-nearest-neighbor interaction and the in- 



terchain interaction. The former is indicated by powder 
neutron scatteringi^ and ESR — measurements and de- 
cisively confirmed by the consideration of Schottky heat 
capacity. ^'^ The latter is confirmed by careful analysis of 
powder neutron scattering measurements »ii In the anti- 
ferromagnetic long-range-ordered state, spins adjacent to 
impurity must be parallel as is schematized in Fig.n(b). 
Therefore, ferromagnetic next-nearest-neighbor interac- 
tion favors antiferromagnetic ordering. Interchain inter- 
action is, of course, necessary for the three-dimensional 
ordering due to Mermin- Wagner theorem,-^ 

Since the low energy excitation in doped NENP was 
well described by the effective Hamiltonian Eq. a 
similar approach could be useful to understand the or- 
dered state of doped PbNi2V2 08. Effective Hamiltonian 
for Mg2+ {S = 0)-doped PbNi2V2 08 will be 

Hi — Jnnnsi • S2, Jnnn < 0, (3) 

which is actually the Hamiltonian for ferromagnetic 
dimer. The ground state is 5 = 1 triplet states, where 
the spins are parallel to each other, and favors three- 
dimensional ordering. On the other hand the effective 
Hamiltonian for PbNi2V208 doped with a magnetic im- 
purity is expressed as three spin ring as is shown in Fig.^ 

(c); 

^.2 = Jai{si ■ Sm + Sm ■ •S2) + JnNNSi ■ S2. (4) 

In this case the ground state depends on the sign 
of Jm and also the type of spin interaction such as 
Heisenberg, Ising, and XY. Because they are subject 
to the species of impurity ions M, systematic study on 
doped PbNi2V2 08 may be valuable to reveal the nature 
of the interaction among the impurity spins and edge 
spins. Although there have been several studies in non- 
magnetic impurity-doped PbNi2V208r2*i2*i^'i^*i^ study 
of magnetic- impurity doped PbNi2V2 08 is very rare and 
even preliminary.^^ We will report a complete study in 
Pb(Nii_^M:r)2V208 (M = Mn, Co, Cu, and Mg) and 
have found the drastic dependence on the kind of im- 
purities in the temperature - concentration (T-x) phase 
diagram. 

To study the magnetic excitations neutron inelastic 
scattering would be the most effective tool if a sin- 
gle crystal of reasonable size were available. However, 
PbNi2V2 08 is incongruent and it is hard to obtain large 
single crystal. Therefore the best we can do is bulk 
measurements such as magnetic susceptibility and heat 
capacity on polycrystalline samples. Experimental de- 
tails will be briefly mentioned in section II. In section 
HI A and C the value and number of effective spins in- 
duced by impurities are carefully discussed and the sign 
of Jm are obtained. The value of effective spin is es- 
timated to be \Sai - 1| in Pb(Nii_:rMj,)2V208 (M = 
Mn, Co, and Cu) and Jm is proved to be antiferro- 
magnetic. In section HI B antiferromagntic transition 
is confirmed in Pb(Nii_a;Cua;)2V208 at extremely low 
temperatures. In section HI D the T-x phase diagrams 



are obtained in Pb(Nii_::cM^)2V208 (M = Mn, Co, Mg, 
and Cu). In section IV the weak ferromagnetism in 
Pb(Nii_j:Coj:)2V208, the drastic dependence of the tran- 
sition temperatures on the species of the impurities in 
T-x phase diagram, and the comparison with impurity- 
doped spin-Peierls material CuGeOa are discussed. 

II. EXPERIMENTAL DETAILS 

Polycrystalline samples of Pb(Nii_2:A^a;)2V208 {M — 
Mn, Co, Mg, and Cu) were prepared by solid state reac- 
tion method. Aligned samples were prepared by apply- 
ing magnetic field to the mixture of powder sample and 
epoxy resin (Stycast 1266) at room temperature^^" DC 
magnetic susceptibilities were measured down to 2.0 K 
in = 0.1 T magnetic field by SQUID magnetometer 
(MPMS-XL of Quantum Design Corp.) AC susceptibili- 
ties were measured down to 50 mK in zero static magnetic 
field by using an adiabatic demagnetization refrigerator 
(/iFridge of Cambridge Magnetic Refrigeration Corp.). 
Heat capacities were measured down to 0.4 K by the re- 
laxation method (PPMS with Helium-3 option of Quan- 
tum Design Corp.) in a magnetic field up to 12 T. 

III. EXPERIMENTAL RESULTS 

A. DC Magnetic Susceptibility 

Figure [3 shows that a magnetic anomaly is detected in 
the DC magnetic susceptibility in Pb(Nii_2,Mn2:)2V208 
(x = 0.050) aligned sample. A sharp peak in || c 
axis and a slight anomaly in ± c at 2.3 K in the 
inset of Fig. [21 confirms that the ground state is an easy- 
axis type antiferromagnetic long-range order. The phase 
transition is observed also in the polycrystalline sample 
of a; = 0.030. No anomaly in the lower concentration 
sample x = 0.010 is simply because of our experimental 
limit. 

Figure 13 shows the temperature dependence of the 
magnetic susceptibilities in Pb(Nii_j:Co2:)2V208. They 
were measured both in a field cooling (EC) and in a 
zero-field cooling (ZFC) processes. The impurity-induced 
anomalies are observed at T = 2.2 K for x ~ 0.010 and 
T = 4.2 K for x = 0.020, and the hysteresis between FC 
and ZFC were detected below the anomaly's tempera- 
ture. The magnetic susceptibilities seem to have a finite 
values for both processes at T = K, which would be 
attributed to an antiferromagnetic transition with weak 
spontaneous magnetization. In mixed spin materials spin 
glass is another possibility but the jump in the heat ca- 
pacity, which we will see later, prefers weak ferromag- 
netism. The weak ferromagnetism is possible owing to 
the existence of the Dzyaloshinskii-Moriya antisymmet- 
ric interactioniSiSS as will be discussed in the section Hvl 

No anomalies are detected in any concentration sam- 
ples of Pb(Nii_a;Cua;)2V208 in the measurements down 
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FIG. 2: Temperature dependence of the magnetic susceptibil- 
ities in Pb(Nii_2,Mn2,)2V208. Inset shows the anisotropy of 
the aligned poly crystalline sample with x = 0.050. 



FIG. 4: Temperature dependence of the magnetic susceptibil- 
ities in Pb(Nii_3;Cu^)2V208. 
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FIG. 3: Temperature dependence of the magnetic susceptibili- 
ties in Pb(Nii_a:Coi:)2V208. The inset is the low temperature 
susceptibility in a; = 0.01 sample. 



to T = 2 K. In contrast with other impurities Cu-doping 
only enhances Curie-like susceptibility at low tempera- 
tures which is shown in Fig. 0] 

Hereafter we will estimate the number and the value 
of impurity-induced effective spins based on VBS model, 
where the magnetization in a finite chain is presumed 
to be induced only by the edge spinsn^^ Let us assume 
that the susceptibility is expressed as a sum of the two 
contributions in the low impurity concentration region: 

(5) 

where Xpuic is the susceptibility in pure PbNi2V2 08 and 
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FIG. 5: Curie constants of Pb(Nii_i,Cui)2V208 from the 
measurement of the magnetic susceptibility and the calcu- 
lated values from the three models. 



Xpara is the additional susceptibility due to impurity- 
induced spins. If Xpara obcys Curie- Weiss law, the num- 
ber and the value of the impurity-induced-effective spins 
will be obtained from Curie constants. As was already 
mentioned in the introduction the ferromagnetic next- 
nearest-neighbor interaction exists, which is in fact con- 
firmed by the observation of effective spin with = 1 in 
Pb(Nii_j;Mgj;)2V208 by some of the present authorspi^ 
Therefore, the following three models are possible in 
impurity-doped PbNi2V2 08; 

(i) The interaction between the edge spin and the im- 
purity spin, Jm in Eq. (0}, is very small and edge 
spins are coupled by ferromagnetic next-nearest- 
neighbor interaction. One effective spin with S = 
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FIG. 6: Curie constants of Pb(Nii_^Mn^)2V208 from the 
measurement of the magnetic susceptibihty and the calculated 
values from the three models. 

1 and one with S = Sm will be detected. 

(ii) Jm is ferromagnetic and one effective spin with S = 
Sm + 1 will be detected. 

(iii) Jm is antiferromagnetic and one effective spin with 
S = \Sm — M will be detected. 

Curie constants of Pb(Nii_a;Cua;)2V208 are plotted as 
a function of x in Fig.O Three lines are calculated values 
from the above three models. The experimental data are 
close to the theoretical line of the model (iii) and this 
suggests that Jcu in Eq. 10} is antiferromagnetic. 

In the same way Xpara of Pb(Nii_a;Mna;)2V208 were 
fitted by the Curie- Weiss law and the experimental 
Curie constants and the calculated values of the three 
models are plotted in Fig. |^ The same argument on 
Pb(Nii_2;Cua;)2V208 may lead to the conclusion that 
the most suitable theoretical value is the model (iii), at 
least for a; — > 0. Hence the impurity-induced spins have 
S = 3/2 degree of freedom at each impurity site and Jun 
is antiferromagnetic. 

We assumed g — 2.2 implicitly because Zorko 
et alm^ has shown that the broad ESR resonance has 
g ~ 2.2 between room temperature and 70 K for 
Pb(Nii_3;Mg3;)2V208 {x < 0.24). However, we have to 
be careful about it because the ESR resonance was at- 
tributed to the bulk spins (on Ni^+ ions) and not that 
on impurity ions. At low temperatures they observed the 
shift of the resonance line, which can be due to the grow- 
ing contribution of impurity-related resonance There- 
fore it might be difficult to definitely determine which 
of the three models is the best description of the spin 
state only by the susceptibility measurement. Since mag- 
netic entropy also gives information on the value and the 
number of spins and it is independent on g value, heat- 
capacity measurements might be more useful to discuss 



the model precisely. We will see that our argument here 
is consistent with heat capacity measurements in III C. 



B. AC Susceptibility at Very Low Temperatures 

Although the DC magnetic susceptibility and heat ca- 
pacity in Pb(Nii_2;Cua;)2V208 have been measured down 
to 2 K in the same way as the samples doped with other 
species of impurities, three-dimensional ordering has not 
been found ye^ (see also lIII A|) . There arizes a question 
if there is no transition in Pb(Nii_a;Cu2;)2V208 or the 
transition temperature is extremely low. To search for 
the magnetic phase transition in Pb(Nii_a;Cua;)2V208, 
measurements of AC susceptibility at low temperatures 
(down to 50 mK) were performed using the adiabatic 
demagnetization method. The real parts of the AC sus- 
ceptibilities in Pb(Nii_a;Cua;)2V208 at very low tempera- 
tures are shown in Fig.[7| here the absolute value was not 
exactly obtained because of the specification of the AC 
susceptibility method but it is sufficient to determine the 
existence or absence of the transition and also to deter- 
mine the transition temperature if it exists. The external 
static magnetic field was not applied in this experiment. 

For the sample with x = 0.010, the change of the sus- 
ceptibility is likely to follow the Curie- Weiss law and 
no anomaly was found above 50 mK. For the sample 
with X = 0.020, however, an anomaly appears in the 
susceptibility-temperature curve at ^ 0.2 K and the 
susceptibility-temperature curve becomes flat below this 
temperature. For the samples with x > 0.30 we observed 
a distinct peak in the susceptibility-temperature curve 
and the peak temperature increases with the impurity 
(Cu^"*") concentration. 

This behavior is qualitatively the same as that 
of Pb(Nii_a;Mgj;)2V208, which has three-dimensional 
ordering at low temperatures, and we attributed 
the peaks in Pb(Nii_a;Cua;)2V208 to magnetic order- 
ings. Hence, we concluded that the ground state of 
Pb(Nii_j;Cua;)2V208 is also the ordered state as the sam- 
ples doped with other impurities. However, temperature 
scale of the ordering temperature is quite different. This 
will be discussed in section IV. 



C. Heat Capacity 

We obtained the magnetic heat capacity Cm 
by subtracting the heat capacity of PbMg2V208, 
which is isostructural to PbNi2V208 from that of 
Pb(Nii_a;Ma;)2V208. The magnetic entropy was ob- 
tained according to the equation S^i = Jq ^dT. Figure 
IHl shows the magnetic heat capacity (upper panel) and 
magnetic entropy (lower panel) of Pb(Nii_j;Cua;)2V208 
with X = 0.050. For _ff = T a small anomaly is ob- 
served at ~ 0.5 K, which is probably attributed to the 
antiferromagnetic ordering because the temperature is 
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FIG. 7: Temperature dependence of the AC susceptibilities 
in Pb(Nii_a;Cua;)2V208 at very low temperatures. 

close to the peak temperature, 0.46 K, in the AC sus- 
ceptibility. This anomaly diminishes by the application 
of a magnetic field. This is the typical behavior of classi- 
cal antiferromagnetic transition. In magnetic field there 
appears a Schottky-like broad anomaly at higher temper- 
ature and the anomaly's temperature increases with the 
magnetic field. This behavior was already reported in 
Pb(Nii_j;Mgj;)2V208 (Ref.llll) and the difference is the 
energy scale of the ordered state such as the transition 
temperature or the critical field for the disappearance of 
the ordered state. 

If we assume that the system can be approximated by 
an assembly of the free spins (effective spins) of the same 
S value, the entropy of the system will be saturated at 
Sni — Nks \n{2S + 1). Note that the above expression 
does not contain the g value. Therefore the magnetic en- 
tropy would be more appropriate quantity to obtain the 
value and the number of spins than the Curie constant. 
The saturation values of the three models which is de- 
scribed in section III A are shown in the lower panel of 
Fig. |Hl(solid fines). We can see that the entropies at the 
various values of the magnetic fields converge to almost 
the same values at T 12 K and a point of inflection 
seems to exist there. 

In addition to the induced moments as effective spins, 
there exist other contributions to magnetic entropy. 

(a) One is the contribution from the spin-gap ex- 
citation. Haldane-like gap excitation, in which 
the energy scale is about 20 K, was observed in 
Pb(Nii_3:Mg3:)2V208 by the high field magneti- 
zation measuremeniiSi and coexistence of gapless 
and gap excitations, which was originally discov- 
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FIG. 8: Magnetic heat capacity (upper panel) and magnetic 
entropy (lower panel) of Pb(Nii_iCua;)2V208 with x = 0.050. 

ered in doped spin-Peierls material CuGeOa, was 
suggested in the ordered phase. Although the con- 
tribution from the Haldane-like gap excitation is 
small at low temperatures, it increases with tem- 
perature. 

(b) Another possible source of the error may come 
from the subtraction of lattice heat capacity; the 
assumption that heat capacity in PbMg2V208 is 
lattice heat capacity in Pb(Nii_a:-^^a;)2V208 might 
not be perfect. 

Compared with the contribution of (a), that of (b) would 
be negligible in the temperature regions of Figs. |H1 be- 
cause the lattice heat capacity changes as ~ T"^. Hence 
we can simply assume that the magnetic entropy is the 
sum of contribution from the effective spin and spin-gap 
excitations, where the former is dominant in low tem- 
perature range and the latter is dominant in high tem- 
perature range. In Fig. [S] entropies with various Cu^+ 
concentration seem to converge near 12 K and simulta- 
neously the entropies with low Cu^"*" concentration have 
a point of inflection near the same temperature. Consid- 
ering that the spin gap in PbNi2V208 is about 20 K, the 
above experimental facts suggest the saturation of the 
entropy due to impurity-induced effective spins and the 
onset of contribution from the spin-gap excitation near 
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FIG. 9: Magnetic heat capacity (upper panel) and mag- 
netic entropy (lower panel) of Pb(Nii_2,Mna;)2V208 with 
X — 0.010. The inset shows AC susceptibility in the same 
sample measured in the adiabatic demagnetization refrigera- 
tor. 



12 K. The experimental value of the entropy at T 12 
K is about 0.7 J Pb-mol~"'^K~^ and this is the closest to 
the saturation entropy, iVfcb ln(2 x i + 1) = 0.576 J Pb- 
mol"^, calculated based on model (iii) proposed in the 
section III A. Hence, we conclude that two edge spins and 
the spin of impurity Cu^^ ion are antiferromagnetically 
coupled and one effective spin of 5 = 1/2 appears near 
the impurity spin. The conclusion is consistent with the 
argument in the section III A. 

Magnetic heat capacities and magnetic entropies 
of Pb(Nii_:rMn:r)2V208 with x = 0.010 and 
Pb(Nii_^Co^)2V208 with x = 0.010 are shown in Figs.EI 
and 1101 respectively. For Mn-doped sample, an anomaly 
due to the antiferromagnetic ordering is found in _ff = T 
at 1.1 K, which is the same temperature of the anomaly 
in the AC susceptibility (see the inset). The anomaly 
disappeared above the magnetic field oi H = 2.0 T. For 
Pb(Nii_2,Coa,)2V208 a broad peak of the magnetic heat 
capacity in zero magnetic field is observed at T = 2.1 K, 
which is the same temperature of the anomaly in the 
magnetic susceptibility (see Fig. O and the inset.) The 
weak ferromagnetism is concluded from both heat capac- 



FIG. 10: Magnetic heat capacity (upper panel) and magnetic 
entropy (lower panel) of Pb(Nii_i,Co2,)2V208 with x — 0.010. 



ity and the magnetic susceptibility measurements. The 
anomaly is suppressed by the application of a magnetic 
field and disappears at H — 4.0 T. 

Magnetic heat capacity and entropy in both 
Pb(Nii_,Mn^)2V208 and Pb(Nii„,Co,)2V208 are 
qualitatively the same as in Pb(Nii_2;Cua;)2V208 and the 
similar argument leads to the same conclusion; the model 
(iii) proposed in the section III A is the best to describe 
the effective spin induced by impurities. Jm in the ef- 
fective Hamiltonian in Eq. Q is antiferromagnetic and 
one effective spin with = jS'jv/ — 1| is detected in the 
magnetic-impurity-doped PbNi2 V2 Os . 



D. T—x phase diagram 

T-x phase diagrams are qualitatively independent on 
the species of impurities as we can see in Fig. 1111 . The 
transition temperatures increase with the impurity con- 
centration and they reach their maximum at some con- 
centrations. In higher concentration region the mod- 
est decreases are observed on Pb(Nii_j:Ma;)2V208 {M = 
Mg and Co) where the higher concentration samples are 
available. We find the similar feature of the T-x phase 
diagram in other spin-gap materials such as spin-Peierls 
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CuGe03f2ii22i2^ two-leg spin ladder SrCusOgf^i and in- 
teracting dimer TlCuClsi^ This shows that the phe- 
nomenon, i.e., impurity-induced antiferromagnetic long- 
range order in spin-gap materials, is universal as already 
discussed briefly in Refs. .ICk .26b. i27. and .28i. 

We find another common feature in the lightly-doped 
region. Small amount of impurities induce antiferromag- 
netic long-range order and the transition temperature 
does not decrease drastically in log-log scale, which is 
shown in Fig. ^] (b), and this suggests no threshold 
concentration in Pb(Nii_j;Afa;)2V208. The absence of 
the threshold was also reported in Zn-doped CuGeOa by 
Manabe et alm^ They studied the T-x phase diagram in 
low-concentration region and showed that the relation 



Tjsi = A exp ( 

X 



(6) 



holds in very low-temperature (T < 8 x 0.8 K) and low- 
concentration (a; < 5 X 10~^) region. 

The unique feature in Pb(Nii_3;Mj.)2V20g is that the 
transition temperatures drastically depend on the species 
of the impurities. The maximum transition temperature 
of Pb(Nii_a,C0a;)2V208 is more than 10 K, while that of 
Pb(Nii_a;Cua;)2V208 is less than 1 K. Non-magnetic im- 
purity Mg^"*" induces lower transition temperature than 
Co^^ but, somehow, higher temperature than Mn^+ or 
Cu^+. The seemingly mysterious behavior can be ex- 
plained by considering effective Hamiltonian in Eqs. (jS)) 
and Q . The impurity dependence of the transition tem- 
peratures will be semiquantitatively discussed in the fol- 
lowing section. 
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IV. DISCUSSION 

The low-temperature phase of Pb(Nii_2:Co2:)2V208 
showed a hysteresis in the susceptibility vs temperature 
curve and it was not a simple antiferromagnetic phase. 
The low-temperature phase of Pb(Nii_a;Coa;)2V208 can 
be attributed to the weak-ferromagnetic phase. The rea- 
son for the occurrence of the weak-ferromagnetic phase 
is due to the structure of PbNi2V2 08, where the Ni 
chain (spin chain) is constructed as the four-fold screw 
chaini^ The screw chain does not have an inversion cen- 
ter between the neighboring Ni sites and therefore the 
Dzyaloshinskii-Moriya (DM) interaction^^'^" exists be- 
tween the two neighboring spins. As is well known, DM 
interaction can cause canting of the sublattice magnetiza- 
tion of the antiferromagnetic phase and the state becomes 
weak-ferromagnetic phase. 

In Fig. ^2 we see that the transition temperatures of 
Pb(Nii_j;Ma;)2V208 differ much depending on the kind 
of impurities M even though the qualitative behavior is 
independent. The transition temperature is obtained as 
the divergence of the staggered susceptibility but there 
is no theoretical calculation for our specific case. How- 
ever, it is possible to semiquantitatively explain the T-x 



FIG. 11: (a) Temperature vs impurity concentration (T- 
x) phase diagram in impurity- doped Haldane material 
PbNi2V208. (b) T-x phase diagram in the lightly doped 
region, where the data are plotted in log-log scales. 



phase diagram such as extremely low transition temper- 
ature in Pb(Nii_a:Cua:)2V208 and very high transition 
temperature in Pb(Nii_2;Co2;)2V208. 

As was already mentioned, the low energy excitation 
in the disordered state in a doped-Haldane chain is well 
explained by the effective Hamiltonian in Eq. jSJl. For 
the appearance of the ordered state in coupled-Haldane 
chains, on the other hand, the interchain interaction must 
be considered in addition. However, the consideration 
of the effective Hamiltonian Eq. ||2Jl which includes only 
Jnnn and Jm would still be meaningful if we assume 
that impurity doping does not affect the interchain in- 
teraction. We could discuss the transition temperatures 
in PbNi2V2 08 doped with various impurities by calcu- 
lating the ground state energy in effective Hamiltonian 
Eqs. Q and and also by calculating the energy in 
local collinear spin structure; parallel in Fig. ^ (b) or 
antiparallel in Fig. ^ (c). 
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Firstly let us examine why Pb(Nii_a;Cua;)2V208 
or Pb(Nii_a;Mn3;)2V208 has lower transition temper- 
atures than Pb(Nii_a;Mg3;)2V208- Effective Hamilto- 
nian for Pb(Nii_2;Mg2;)2V208 is Eq. ^ and the ground 
state is ferromagnetic triplet state, as is visualized in 
Fig. n^b). Such a coUinear state conserves the co- 
herence of antiferromagnetic correlation in the chain 
around the impurities. Hence three-dimensionally or- 
dered state is favored by the local inchain interac- 
tion, ferromagnetic Jnnn, as well as interchain inter- 
action. On the other hand the effective Hamiltonian 
for Pb(Nii_,Cu,)2V208 or Pb(Nii_,Mn,)2V208 is ex- 
pressed by Eq. Q). The local antiparallel spin struc- 
ture (see Fig. n(c)), which favors coherence of antiferro- 
magnetic correlation in the chain, is neither the ground 
state nor an eigenstate of the effective Hamiltonian. 
This means that the three-dimensionally ordered state in 
Pb(Nii_:rCu:r)2V208 or Pb(Nii_xMn:r)2V208 is favored 
by interchain interaction but it is disturbed by the local 
inchain interaction, Jm- Therefore Pb(Nii„^Cua;)2V208 
or Pb(Nii_2;Mn2;)2V208 has lower transition tempera- 
tures than Pb(Nii_^Mg^)2V208. 

Secondly let us discuss the transition temperatures 
in Pb(Nii_,Cu,)2V208 and Pb(Nii_,Mn,)2V208. The 
ground states of both of the effective Hamiltonians are 
not a coUinear spin state. This fact does not favor the 
antiferromagnetic coherence in the isolated chain and the 
appearance of ordered state is the result of the reduction 
in the energy of the ordered state due to the interchain in- 
teraction. In this case the transition temperature may be 
mainly determined by the energy difference between the 
ground state of the effective Hamiltonian and coUinear 
spin state. We define 

A „ _ /-E'collinear ~ Eq\ 

Aii7\/=Cu.Mn = 7 t; , (I) 

^ JaiSIJM /M=Cu,Mn 

where 

-BcoUinoar.Af = (colUnear, M|'H2 |coUinear, M) 
IcolUnear, M) = \sl ^ 1/2, S^i ^ -Sm , 4 = 1/2). (8) 

For the appearance of three-dimensionally ordered state, 
large AEm requires large energy reduction in the or- 
dered state and this results in the lower transition tem- 
perature. The maximum transition temperatures T^'^^, 
Eg,Mi -£'coiiinear,Mj and IS.Em are summarized in TableQ] 
We find Ai?cu > Ai^Mn and this is consistent with lower 
transition temperature in Pb(Nii_2;Cua;)2V208 than in 
Pb(Nii_,Mn,)2V208. 

Thirdly let us discuss the highest transition tem- 
peratures in Pb(Nii_a;Coa;)2V208. Since we have 
already explained that Pb(Nii_2;Mg2;)2V208 has the 
highest transition temperature among Mg-, Cu- 
, and Mn-dopcd PbNi2V208, the comparison be- 
tween Pb{Nii_2,Co2,)2V208 and Pb(Nii_:rMga;)2V208 is 
enough. So far we assumed implicitly that the spin in- 
teractions are isotropic. In some magnetic ions such as 
Co^"'", however, the interaction is possibly Ising-like be- 
cause of remnant orbital momentum. In this case the 



first and the second term in Eq. @ are modified; 

Ha = Jco{sl ■ S'co + S'co ■ si) + JnnnSi • S2- (9) 

The ground state of the effective Hamiltonian is spin 
coUinear structure, as is similar to Pb(Nii_2;Mg3;)2V208, 
and the lower the ground state energy is, the higher the 
transition temperature will be. Useful physical quanti- 
ties for Pb(Nii_,Mg,)2V208 and Pb(Nii_,Co,)2V208 
are summarized in Table ^ The ground state en- 
ergy in Pb(Nii_a;Co2:)2V208 is lower than that of 
Pb(Nii_^Mg^)2V208 by -(3/2)Jco and, therefore, the 
former has higher transition temperature. This means 
that Pb(Nii_j;Co2;)2V208 has the highest transition tem- 
perature. 

The T-x phase diagram in Pb(Nii_j;Ma;)2V208 has 
been semiquantitatively explained by the simple effec- 
tive Hamiltonians based on VBS model. The drastic 
difference between magnetic impurity (Cu^+ or Mn^+) 
and non-magnetic impurity (Mg^+) is due to the differ- 
ent ground states. The highest transition temperature 
in Pb(Nii_a;Co2;)2V208 could be ascribed to Ising-type 
interaction. The effective Hamiltonians were successfully 
applied to Haldane materials but they can not be applied 
to other spin-gap materials such as S = 1/2 spin-Peierls 
CuGeOa. Indeed there is no drastic difference in the 
transition temperatures between non-magnetic impurity 
(Mg^+ or Zn^+) and magnetic (Ni^+ 5* = 1) impurity- 
doped CuGeOs.^^ 

In doped CuGeOs some of the present authors found 
that there is a first-order phase transition between 
the two kinds of antiferromagnetic phases with the 
change of the concentration of non-magnetic impu- 
rity Mg2+ (Refs. ^ and ^ or Zn^+iS^ Two phases 
are the dimerized- antiferromagnetic (DAF) phase and 
uniform-antiferromagnetic (UAF) phase. This phe- 
nomenon has been studied in detail by the various meth- 
ods; thermal conduct ivityjS^ Raman scatteringf^^ neu- 
tron diffraction,^^^^ pressure effect,^^'^^^ ESR,^ and 
synchrotron x-ray scattering4i In Fig. Illf al we see that 
the antiferromagnetic phase of Pb(Nii_a;Mga;)2V208 ex- 
tends very widely to x ^ 1 (PbMg2V208 is not mag- 
netic) and there seems no compositional phase tran- 
sition in this wide region. On the contrary in Mg- 
doped CuGeOs the concentration region of the antifer- 
romagnetic phase is more limited and there is a first- 
order compositional phase transition as mentioned above. 
The concept of the "impurity-induced antiferromagnetic 
phase" may only be applied to the lightly doped region 
in Pb(Nii_a;Mga;)2V208. When the concentration x of 
Mg is very large, magnetic Ni^+ ions are randomly dis- 
tributed three dimensionally among the non-magnetic 
Mg^+ ions and the ordered state may be constructed 
among these diluted Ni^+ spins. In CuGeOa the spin- 
gap state is caused by the spin-lattice interaction. In 
CuGeOa and doped CuGeOs spin gap is created by the 
spin-Peierls transition and the antiferromagnetic phase in 
the low impurity concentration region (DAF phase) has 
two order parameters: dimerization and antiferromag- 
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TABLE I: Summary of the maximum transition temperature, the ground state energy, colhnear spin state energy, and the 
scaled energy difference. Details are explained in the text. 

M = Cu M = Mn M = Mg M = Co 

IW^ 0.65 K{x = 0.1) 2.2 K (a; = 0.05) 3.5 K (a; = 0.08) 10.4 K (a; = 0.15) 

n_/ -u ■xj ni nr 

rt /t2 tt'z rti rtj, 

Eg,m —Jcu + 1/4Jnnn — 7/2JMn + 1/4Jnnn 1/4Jnnn — 3/2Jco + 1/4Jnnn 

-Bcollinear.M — 1/2 Jcu + 1/4 JnNN — 5/2JMn + 1/4Jnnn 1/4Jnnn — 3/2Jco + 1/4Jnnn 

AEm 2 4/5 



netic long-range order. The difference of DAF and UAF 
phases in doped CuGeOs is the existence or absence of 
the lattice distortion (dimerization) . On the contrary in 
PbNi2V2 08 the spin gap is intrinsic and in the impurity- 
induced antiferromagnetic phase there is ony one order 
parameter (antiferromagnetic long-range order) . The dif- 
ference of the numbers of the order parameters may cause 
the different behavior of the compositional change be- 
tween Cui-^Mg^GeOs and Pb(Nii_^Mg^)2V208. 

V. SUMMARY 

Impurity-induced three-dimensional ordering is 
systematically studied in doped Haldane material 
PbNi2V208 by use of DC and AC magnetic suscepti- 
bility, and heat capacity, on many samples with various 



species of impurities and concentrations. Complete T-x 
phase diagrams are obtained for each impurity and 
qualitatively common features to the doped spin-gap 
materials are observed. The unique feature is found in 
the drastic dependence of the transition temperatures 
on the species of the impurity, which is explained by 
effective Hamiltonian based on VBS model. 
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